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Using the 2dF galaxy redshift survey to detect 
gravitationally-lensed quasars 
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ABSTRACT 

Galaxy redshift surveys can be used to detect gravitationally-lensed quasars if the 
spectra obtained are searched for the quasars' emission lines. Previous investigations of 
this possibility have used simple models to show that the 2 degree Field (2dF) redshift 
survey could yield several tens of new lenses, and that the larger Sloan Digital Sky 
Survey should contain an order of magnitude more. However the particular selection 
effects of the samples were not included in these calculations, limiting the robustness 
of the predictions; thus a more detailed simulation of the 2dF survey was undertaken 
here. The use of an isophotal magnitude limit reduces both the depth of the sample and 
the expected number of lenses, but more important is the Automatic Plate Measuring 
survey's star-galaxy separation algorithm, used to generate the 2dF input catalogue. It 
is found that most quasar lenses are classed as merged stars, with only the few lenses 
with low-redshift deflectors likely to be classified as galaxies. Explicit inclusion of these 
selection effects implies that the 2dF survey should contain 10 lenses on average. The 
largest remaining uncertainty is the lack of knowledge of the ease with which any 
underlying quasars can be extracted from the survey spectra. 

Key words: gravitational lensing - galaxies: surveys - methods: data analysis. 



1 INTRODUCTION 

Gravitationally-lensed quasars are extremely valuable, and 
so a great deal of effort is expended searching for new lens 
systems. Most lens surveys are based on high-resolution re- 
observation of known quasars, but lensed quasars can also 
be found galaxy redshift surveys (GRSs) if the spectra ob- 
tained are searched for the presence of quasar emission lines. 
As well as being efficient, this form of lens search is comple- 
mentary to the conventional surveys, being primarily sensi- 
tive to low-redshift deflectors - the galaxies that would enter 
the GRS a priori. This is important as the proximity of the 
lens galaxy can facillitate a number of unique observations, 
as exemplified by Q 2237+0305 (Huchra et al. 1985), the 
only quasar lens yet found in a GRS. 

A more detailed discussion of these points is contained 
in Mortlock & Webster (2000b), one of the two previous in- 
vestigations of quasar lensing in redshift surveys. An 'initial, 
order of magnitude estimate' of the number of lenses ex- 
pected in future GRSs by Kochanek (1992) suggested there 
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would be one lens per ~ 10 5 galaxy spectra, but this was 
shown to be an under-estimate by up to an order of mag- 
nitude, as the combined light of the quasar images and the 
lens galaxy must be included in the calculation (Mortlock & 
Webster 2000b). Both these these results were obtained for 
an idealised GRS, characterised only by the number of spec- 
tra obtained and a survey's magnitude^ limit. It is important 
to test the robustness of these simple predictions, which can 
be done most profitably by simulating a real survey in de- 
tail. The two obvious candidates are the 2 degree Field (2dF) 
galaxy survey (e.g. Colless 1999; Folkes et al. 1999) and the 
Sloan Digital Sky Survey (SDSS; e.g. Szalay 1998; Loveday 
& Pier 1998). The former was chosen primarily as it subject 



t All magnitudes are in the Bj system, but this subscript is omit- 
ted for brevity. Also, m denotes total apparent m agni tude, m' 
isophotal apparent magnitude (as defined in Section 2.3), and M 
denotes absolute magnitude. An Einstein-de Sitter cosmologcal 
model is assumed for simplicity; the choice cosmology does not 
have a significant impact on the inferred properties of the 2dF 
galaxy survey or the lens statistics, as the galaxies are so nearby 
(Kochanek 1992). Finally, Hubble's constant is taken to be 70 km 
s _1 Mpc -1 . 
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to a number of comlex, surface brightness-related selection 
effects, whereas the SDSS is much closer to the ideal galaxy 
sample modelled in Mortlock & Webster (2000b). Further- 
more, the 2dF survey will be completed first, offering the 
most immediate opportunity to discover new quasar lenses 
in a galaxy survey. 

The selection criteria of the 2dF GRS and its parent sur- 
vey are carefully defined in Section ^. The main aim of these 
selection criteria is to exclude stars (Section ||) whilst max- 
imising the completeness of the galaxy sample (Section^), 
although care was also taken to reject unwanted merged im- 
ages (Section - In order to determine the number of lenses 
in the survey, the same selection effects must be applied to 
a simulated lens population, and the spectral sensitivity of 
the lens search must also be included (Section |^). The re- 
sults of the this calculation and the relative importance of 
the various survey characteristics are given in Section ^. The 
main conclusions are summarised in Section H. 



2 THE 2DF GALAXY REDSHIFT SURVEY 



The 2dF instrument (Section 2.1) is used for a number of 



proj ects . The largest is the eponymous galaxy survey (Sec- 
tion 2.2), the properties of which are, however, determined 



mainly by its parent sample (Section 2.c ) 



2.1 The 2dF instrument 

The 2dF instrument (Taylor 1994; Cannon 1995) is a 
wide-field, multi-fibre spectroscopic survey facility for the 
Anglo- Australian Telescope (AAT). Corrective optics give 
aberration-free images over the full 2 deg-diameter field 
of view, but the main innovation is the automatically- 
configured 400-fibre spectrograph in the prime focus plane 
of the telescope. The positional accuracy of 0.1 arcsec is re- 
quired because the fibres have an effective angular radius 
of only Of — 1 arcsec. This aperture size is comparable to 
the seeing at the site - which varies between 1 arcsec and 
5 arcsec, with a median value of ~ 2 arcsec (Ryan & Wood 
1995) - and so the spectra obtained are inevitably of variable 
quality. 



2.2 The 2dF survey 

The 2dF GRS (e.g. Colless 1999; Folkes et al. 1999) is the 
largest single project that is proposed for the 2dF instru- 
ment]^]. The survey will require ~ 100 nights of observing 
time to obtain the redshifts of iVtot = 2.5 x 10 J objects, 
about 95 per cent of which are expected to be galaxies. The 
spectrograph being used has a resolution of 9 A from ~ 3700 
A to ~ 8100 A, and 45-minute integrations should yield 
spectra with a signal-to-noise ratio of ~ 20 per pixel at the 
survey magnitude limit of ml ~ 19.5. 

In the absence of any other selection effects, the sur- 
vey limit, together with JVtot, would completely characterise 
the survey, and it could then be simulated using standard 

t The 2dF quasar survey (e.g. Boyle et al. 1999a, b) is a sep- 
arate project, but data for both it and the GRS are obtained 
simultaneously. 



techniques. However, the catalogue of galaxy candidates is 
taken from the Automatic Plate Measuring (APM) galaxy 
survey, which is subject to a number of complex, surface 
brightness-related selection effects. 

2.3 The APM galaxy survey 

The APM system (e.g. Kibblewhite et al. 1984), combines 
a high-speed microdensitometer with a dedicated computer 
system capable of on-line image detection. Used to analyse 
United Kingdom Schmidt Telescope (UKST) photographic 
plates, it can provide spatial information on a pixel scale 
of ~ 0.5 arcsec, and relative photometry over several orders 
of magnitude in flux. Scanned at this resolution, a typical 
UKST plate represents ~ 3 GB of information, but the im- 
age detection algorithms allow it to be stored as a ~ 3 MB 
list of image parameters with minimal information loss. The 
APM software defines an image as any connected set of 16 
or more pixels which are all above the surface brightness 
limit. This limit varies between /Ltii m = 24 mag arcsec -2 
and /iun, = 25 mag arcsec -2 . The minimum image area of 
~ 4 arcsec 2 then implies an isophotal magnitude limit of 
ml ~ 22. For each image, 15 parameters (defined in Maddox 
et al. 1990a) are calculated from the sky-subtracted plate 
densities; these comprise the standard output of the APM 
system. 

The APM galaxy survey (Maddox et al. 1990a; Maddox, 
Efsathiou & Sutherland 1990b; Dalton et al. 1994; Maddox, 
Efsathiou & Sutherland 1996; Dalton et al. 1997) is a uni- 
form sample of ~ 2 x 10 6 galaxies to a magnitude limit of 
ml — 20.5. It was compiled from APM scans of 185 contigu- 
ous UKST survey plates, and the galaxy sample was thus 
selected from ~ 2 x 10 7 images brighter than ml ~ 22 spread 
over a total area of 4300 deg 2 . 

2.3.1 Star-galaxy separation 

Given that only a third of objects brighter than m ~ 20.5 
are galaxies, one of the most important aspects of the APM 
galaxy survey is the star-galaxy separation technique. In 
Mortlock & Webster (2000b), image ellipticity was used as a 
generic method of discriminating between point-like images 
and extended sources - this was useful as a simple estimate, 
but more sophisticated methods are required to construct 
a galaxy catalogue from real two-dimensional data. Most 
star-galaxy separation techniques (e.g. Jarvis & Tyson 1981; 
Godwin, Metcalfe & Peach 1983; Maddox et al. 1990a) are 
based on two facts: a galaxy has a lower peak surface bright- 
ness than a star of the same magnitude; and the surface 
brightness profile of a star is determined purely by obser- 
vational effects. Stars and galaxies inhabit different regions 
of the phase-space defined by magnitude, area, peak surface 
brightness, etc., and scatter plots of the data can be used to 
separate the two populations. 

The APM galaxy sample was defined in a more complex 
manner, using the areal profiles of each object as the main 
diagnostic. Amongst the image parameters provided by the 
APM is the area of each image above eight plate density 
levels, Ai, where 1 < i < 8. The plate density is propor- 
tional to (logarithmic) surface brightness fainter than ~ 20 
mag arcsec -2 (Maddox et al. 1990a), so the levels can be 
approximated by 
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p-% = Hum + i A/i, 



(1) 



where /in 1 1 1 

-2 



24.5 mag arcsec - and A/i ~ 0.75 mag 
arcsec ". This approximation slightly underestimates us, 
due to saturation effects, and also ignores plate-to-plate vari- 
ations. The random errors in the AiS were found to be well 
approximated by the Poisson error in the number of pixels 
above the relevant surface brightness. Given the APM scale 
of 0.5 arcsec per pixel, AAi ~ 2\fA~i if Ai is expressed in 
arcsec 2 . The areal profile of an image can be thought of as 
a plot of area versus surface brightness, which is similar to 
a conventional surface brightness profile. As shown in Mad- 
dox et al. (1990a), stars have a characteristic areal profile at 
any given magnitude, which is distinct from the more varied 
profiles of galaxies. However this differentiation breaks down 
for m' > 20 (as there is limited information for such faint 
images) and for m! Ss 17 (as the stars appear extended due 
to saturation effects). 

As a source of extra information, the peak surface 
brightness of each image, /i pca k, and the radius of gyration, 
#G, were also used, giving a total of 10 parameters for the 
star-galaxy separation. The radius of gyration is defined as 



(2) 



where the expectation values are given by surface brightness- 
weighted integrals over the area of the image for which /i < 
jUii m . The uncertainties in /i poa k and 8g are not given in 
Maddox et al. (1990a), but are estimated to be ~ 10 per 
cent. As with the AiS, a distinct stellar locus is apparent for 
both /ipeak and 8q as a function of magnitude - the stars 
have higher peak surface brightnesses and smaller radii of 
gyration than the galaxies. 

The final image classification statistic used by Maddox 
et al. (1990a) combines all 10 image parameters in a x 2 -like 
formulation. For an object of isophotal magnitude m' , it is 
defined as 



V>APM 



2000 



+ 



0c, - 6> G , s (m') 
^200,3 (m') 



i=l 

2 



A itS {m') 



y/lAi^m') 



(3) 



/peak ,/peak,s(^' ) 
\//pcak,s(m') 



where the Ai are in units of arcsec 2 , 8q is in units of arc- 
sec, and /p C ak is the peak surface brightness in linear units, 
defined to match the APM plate density. The 's' subscript 
denotes the parameter value of an isolated star of the same 
isophotal magnitude as the image in question. Each term is 
normalised by the expected uncertainty in the stellar locus, 
as in a \ 2 formulation. Note that the definition of i/'apm 
given in Maddox et al. (1990a) does not include the over- 
all normalisation (i.e. the '1/10' term); it is included here 
to ensure that stellar locus is centred on tpAPM = 0. The 
prefactor is changed from 1000 to 2000 to compensate for 
this. 

Fig. l] shows a plot of i/>apm as a function of isophotal 
magnitude for stars, galaxies and mergers, and is directly 
comparable to Fig. 10 of Maddox et al. (1990a). Some of 
the details differ (e.g. the exact position and width of the 
stellar locus and absence of saturation effects in the simula- 
tion), but the generic features are the same. The stars (See 
Section ^.) are centred on "0apm = 0, and their spread is 
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Figure 1. The primary APM star-galaxy separation parame- 
ter, i/>APM; as a function of APM isophotal magnitude, ml . The 
stars (points) lie along i/'apm = 0, by definition, and the galaxies 
(also points) comprise the well-defined 'ridge' with i/'apm > 1000 
(above the horizontal dashed line), although the distinction be- 
tween the two populations is less clear at fainter magnitudes. Also 
shown are representative populations of lensed quasars (solid cir- 
cles) and pairs of merged images: star-star mergers (open circles); 
star-galaxy mergers (open squares) and galaxy-galaxy mergers 
(open triangles). The vertical dashed line at m' = 19.5 represents 
the 2dF survey limit. 



independent of the point-spread function (PSF) due to the 
normalisations in equation (^). The galaxies (See Section |j|.) 
lie along a 'ridge' with ipAPM ~ 1000, but are less distinct 
from the stars at fainter magnitudes. However ipAPM cannot 
be used to separate merged images from extended objects; 
further specialised statistics are required. 



2.3.2 Merged objects 

Images separated by 5 arcsec are registered as single ob- 
jects by the APM software, and tend to have i/'apm ?Z, 1000, 
as illustrated in Fig. [ll If this was the only star-galaxy sepa- 
ration implemented, most lensed quasars would be included 
in the APM galaxy survey. However Maddox et al. (1990a) 
used the radius of gyration and a saturation parameter to 
detect (and, in general, reject) merged images. 

The radius of gyration, 8g, as defined in equation (^|) is 
a weighted measure of the extent of an image, and is eval- 
uated on a pixel- by-pixel basis by the APM software. If an 
image is assumed to be perfectly elliptical, an approximation 
to the radius of gyration is given by 



J2 7 i=1 fiMAi-A i+1 ) 



1/2 



(4) 



where e is the overall measured ellipticity of the image and 
fi the surface brightness of the ith level in linear units (Mad- 
dox et al. 1990a). For isolated images, 8a — 8' G , with some 
scatter due to random noise; for multiple images #g > 8' G , 
as dumbbell-shaped image pairs have higher second order 
moments than isolated images with the same areal profile. 
From this Maddox et al. (1990a) define 
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Figure 2. A scatter plot of the two parameters used to char- 
acterise merged images in the APM galaxy survey, &apm an d 
MAPM- The dashed line defines the boundary between 'stellar' 
mergers (above the line) and 'galactic' mergers (below the line) 
- only the latter are included in the galaxy survey. The stars 
(points) are centred at &apm = 1 an d /iapm = 1, by definition. 
The galaxies (also points) have &apm — 1 as well, but a broad 
spread in (iapm- Also shown are representative populations of 
lensed quasars (solid circles) and pairs of merged images: star- 
star mergers (open circles); star-galaxy mergers (open squares) 
and galaxy-galaxy mergers (open triangles). 



&APM = 



/o'i 



or 



(5) 



where again the 's' subscript denotes stellar values at the 
same isophotal magnitude as the (merged) object in ques- 
tion. If &apm is much greater than unity the object is quite 
likely to be a merged pair of images or a higher multiple. 

The fc-parameter is sufficient to find pairs of nearby im- 
ages, but is not sensitive to superimposed objects (e.g. a 
star centred on a galaxy) and cannot be used to differenti- 
ate between galaxy-galaxy mergers, star-galaxy mergers and 
star-star mergers. Such a distinction can be made by looking 
at the fraction of the image that is close to saturation on the 
UKST plates. To this end Maddox et al. (1990a) defined 



Capm 



(0.2A 7 + As) /Ai 



[0.2A T , s (m>) + A 8 , s (m')] /Ai, a (m') ' 



(6) 



where the inclusion of A7 (rather than just using A&/A\) 
gives extra stability for fainter images. By definition stars 
(both single and multiple) have /xapm — 1, whereas galaxies 
(and galactic mergers) < /iapm < 1. 

A scatter plot in &apm~/^apm space (Fig. |j| c.f. Fig. 14 of 
Maddox et al. 1990a) can thus be used to separate the var- 
ious types of mergers from isolated stars and galaxies. As 
explained in Maddox et al. (1990a), merged galaxies were 
retained, as were bright galaxies blended with fainter stars, 
but all other stellar mergers were rejected. This was done 
by making a cut 'by eye' in the ^apm-Mapm plane, which is 
shown as the dashed line in Fig. ^. Note also that the nor- 
malisation of both &apm and £*apm is such that there is no 
intrinsic magnitude dependence in the merger classification. 



3 STARS 

As the APM image classification scheme is defined relative 
to the stellar population, stars represent a control sample. 
Misclassified stars are also the primary source of contamina- 
tion in the galaxy catalogue, as they represent such a large 
fraction of sources over the magnitude range of the survey. 

The model of the stellar population used was described 
in Mortlock & Webster (2000b). Briefly, the number counts 
are taken from Bahcall & Soneira (1980) and the images are 
assumed to have a Moffat (1969) PSF, defined by 



/sec (0) — 



(v-m 



i/ v 



i) 



(0 s /2) 2 



(7) 



where 8 S is the full width at half-maximum of the seeing disc 
and 6 is angular position on the sky. As discussed in Mort- 
lock & Webster (2000b), two values of r\ are used: 77 = 2.5, 
the default value assumed by the IMAGE reduction and 
ANALYSIS facility software (Tody 1986); and 77 — > 00, a 
Gaussian PSF. It is important to note that this PSF is for- 
mally defined for an object of a given total, as opposed to 
isophotal magnitude. An iterative procedure is required to 
find the total magnitude, m, which corresponds to a given 
isophotal magnitude m'. The random errors due to the finite 
signal-to-noise of the APM pixels are not explicitly simu- 
lated here; onl y the approximately Gaussian errors discussed 
in Section |2.3| are included. 

Fig. [j]shows the stellar locus as the 'ridge' with t/'apm — 
0. The spread is independent of the PSF as Vapm is nor- 
malised relative to the spread of the stellar images. In re- 
ality the stellar locus is not straight and has non-Gaussian 
outliers (c.f. Fig. 10 of Maddox et al. 1990a), but the most 
important features are present. 

Isolated stars have an even simpler distribution on the 
fcAPM-^APM plot (Fig. ^) used to distinguish between merg- 
ers and isolated objects. Their distribution in &apm is simply 
a Gaussian with spread determined by the assumed error in 



the radius of gyration discussed in Section 2.3. The use of A7 
and A% in the definition of /iapm [equation (6|)] means that 
the fiAPM distribution is marginally magnitude dependent, 
but is also close to a Gaussian for stars with /i pea k — A*8- 



4 GALAXIES 

The main purpose of the morphological classification crite- 
ria of the APM survey is to generate a reasonably complete 
and contamination-free sample of galaxies. The properties 
of the galaxy sample are determined by the intrinsic galaxy 
population (Section WTj) and their observed surface bright- 
ness profiles (Section |4.2| ). The latter dominates both the 
APM classifications (Section fl.3| ) and the 2dF fibre magni- 
tudes (Section 4.4), but the two combine to determine the 
redshift distribution of the survey galaxies (Section 4.5). 



4.1 Populations 

The model of the local galaxy population is described in 
detail in Mortlock & Webster (2000b), and thus only out- 
lined here. Three galaxy types are considered: spirals (S) 
and two classes of elliptical galaxies (E and SO), the relative 
numbers of which are given in Postman & Geller (1984). 
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Each population is assumed to follow the Schechter (1976) 
luminosity functions detailed in Efstathiou, Ellis & Peter- 
son (1988). (Thus M, = -19.4 ± 0.1 for all three types.) In 
the calculation of the galaxies' lensing properties (Sections ^ 
and they are modelled as isothermal spheres (e.g. Turner, 
Ostriker & Gott 1984; Binney & Tremaine 1987), with the 
standard Faber- Jackson (1976) and Tully-Fisher (1977) re- 
lations used to convert luminosities to velocity dispersions. 

4.2 Surface brightness profiles 

In the APM survey, galaxies are distinguished from stars 
primarily by their surface brightness profiles. In particular, 
the ellipticity or elongation of images is not used in the ob- 
ject classification^, and so all galaxies can be assumed to 
have circular isophotes without loss of generality. 

In the absence of atmospheric seeing, the surface bright- 
ness of galaxies is given by their intrinsic (radial) surface 
brightness profiles. Spirals are assumed to follow a Freeman 
(1970) profile, given by 

fa(e) = - e -e/(o.596e g ). (g) 

JS( > 27r(0.596 6» g )2 ' (> 

the two types of ellipticals are modelled by a de Vaucouleurs 
(1948) profile, given by 



= 2P 



M ) = 296.7^ e-W/W 1 '*, 



(9) 



where F is the total flux of the galaxy and 6 = \0\ (where 
6 is the position on the sky relative to the centre of the 
galaxy). Here 9 g — R s /cIa(0, z g ), where <1a is the angular 
diameter distance to the galaxy and R g its intrinsic half-light 
or effective radius. This scales with the galaxy's velocity 
dispersion as 



R& — R 



( a \ Us 



(10) 



where _R g * = (3 ± 1) kpc and u g = 4 ± 1 for Es and SOs 
(Kormendy & Djorgovski 1989) and R g , = (4 ± 1) kpc and 
u g — 3 ± 1 for spirals (de Vaucouleurs & Pence 1978). 

The observed surface brightness profile of a galaxy, 
fg(0), can be found by convolving the intrinsic profile with 
the PSF. Even if the PSF is also rotationally-symmetric, 



(ii) 



/sorB(0')/ S ee ( yj 2 + #' 2 - 299' COS 0J, ) 9'd9' <t> d9' , 

involves a two-dimensional integral, which must be evalu- 
ated numerically. This is time-consuming, and so an approx- 
imation was sought. 

The analytic surface brightness profile used is 

f g (9) = f s (0)eWe-(^^ /e '^ (12) 

£-F 1 ~(V^+^/e' g y 

27<r(2/£) l-P[2/^(9' c /9' g ^f 

where P(z) is the Gamma function and P(z, x) is the incom- 
plete Gamma function. Unfortunately 9' g can only be found 
by solving 



S Ellipticity was the only criterion used for star-galaxy separation 
in Mortlock & Webster (2000b). 



2/e, (yw+oi/e & y 



P [2/e, (9' c /9 g 



1(13) 



numerically, but, more importantly, f g (9) can be inverted, 



0(f) = ^9'i[{9U9' g f + \nU g {Q)/f)] 2n -9' c \ (14) 

provided / < / g (0). The profile can also be integrated, to 
give the flux within a specified angle as 



F e {< 9) = F- 



p 


2/C, (V9*+9?/9' g y 


- P [2/C, {6'je' & f] 


l-P[2/f,(0£/20 g )«] 



(15) 

Each of the three parameters plays a clear role in the 
shape of f s (9). The core width, 9' c , is determined by the 
mainly by the seeing, but has some dependence on the form 
of the PSF. The overall angular scale of the profile, given 
by 9' g , is purely a function of the galaxy's distance and half- 
light radius. For 9 3> 9' c the surface brightness falls off as 
exp[(#/# g )~^], assuming that the faint wings of the galaxy 
are broader than those of the PSFpjj. Hence £ = 1 for spirals 
and ^ — 1/4 for ellipticals, and, if 9' c — 0, equation ( p^[ ) 
reduces to /e and fs. For small 9 equation ( |12[ ) tends to 
/ g oc exp[— (9 2 /29 c 9' g )^], and so the profile has a Gaussian 

core only if £ ~ 1 (i.e., for spirals)]^]. Fig. ^| shows several 
examples of observed galaxies' profiles, comparing the nu- 
merically integrated curves (dashed lines) to the relevant 
form of equation ([H^) for both ellipticals (a) and spirals (b) . 
The agreement for the latter type is excellent, with relative 
errors of no more than 1 per cent for the profiles shown. The 
average error for the elliptical fits is ~ 1 per cent, and peaks 
at ~ 5 per cent for most of the profiles shown. 

Both the approximate depth of the redshift survey and 
the importance of the use of isophotal magnitudes can be 
gauged by plotting the magnitude of a given galaxy as a 
function of redshift (i.e. a Hubble diagram). This is shown in 
Fig. § for both an M» = -19.4 elliptical and an M* = -19.4 



spiral, with the 2dF fibre magnitude (defined in Section 1.4) 
also included. The m(z) curve is the same for the elliptical 
and the spiral as the luminosity evolution of each type is 
assumed to exactly cancel their K- correct ions (See Mort- 
lock & Webster 2000b), but the redshift dependence of the 
isophotal magnitude is somewhat different. For moderate 
redshifts, the isophotal flux of spirals is actually greater than 
that of ellipticals, despite the fact that the latter have more 
peaked surface brightness distributions. This comes about 
as the average flux of the spirals over the region registered 
by the APM is higher. However at high redshifts m' — *• oo 
for all galaxies as they are so faint that their peak surface 
brightness is lower than the APM limit of 24.5 mag arcsec -2 . 



It This is not formally the case for a Moffat (1969) PSF, although 
it only becomes relevant on very large angular scales. 
II An alternative profile with 



fsW = / g (0)e 



(16) 



would have a Gaussian core for all £ with the same properties as 
equation (p^) for large 8, but cannot be integrated analytically. 
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Figure 3. The observed surface brightness profiles of elliptical (a) and spiral (b) galaxies. The dashed lines are the result of convolving de 
Vaucoulcurs (1948) and Freeman (1970) profiles, respectively, with a Moffat (1969) PSF; the solid lines are the analytical approximations 
given in equation (|l^). In both cases the galaxy has absolute magnitude M* = —19.4 and is at a redshift of 0.1, and all the profiles are 
normalised to the peak surface brightness of the intrinsic galaxy profile. From top to bottom the seeing is arcsec, 1 arcsec, 2 arcsec, 3 
arcsec and 4 arcsec. 



This occurs at a higher redshift for ellipticals, as it is deter- 
mined mainly by the central peak of the surface bright ness . 
Over the redshift range of the 2dF survey (See Section 4.5.) 
m' — m ~ 0.5 for most galaxies, and so the redshift coverage 
of the survey is comparable to that of a GRS with a total 
limiting magnitude of mij m ~ 19. 



4.3 APM classifications 

The observed surface brightness profile of a galaxy deter- 
mines its classification as either a stellar or non-stellar 
image. Predictably, almost all the brighter galaxies have 
ipAPM <; 1000, and are classified as non-stellar, as shown 
in Fig. |l|. It is only fainter than the 2dF survey limit of 
m' — 19.5 that the galactic and stellar loci begin to merge. 
The galaxies also occupy a well-defined locus on the ^apm- 
A*apm plot, as shown in Fig. ^| Isolated galaxies are mono- 
lithic and hence have /capm — 0, but the cover a wider range 
in /japm- Fainter and more distant galaxies become progres- 
sively more stellar, and so /japm — > 1, the canonical stellar 
value. Galaxies with ^apm > 0.8 are in fact classified as 
stellar mergers, and rejected from the galaxy sample. Over- 
all, however, it is clear from Figs. |l|and| that very few real 
galaxies with m' < 19.5 are lost from the 2dF catalogue, and 
a Monte Carlo simulation suggests the APM galaxy survey 
has a completeness of C ~ 95 per cent, which is comparable 
to the values estimated by Maddox et al. (1990a). 



4.4 Fibre magnitudes 

In general, the flux from a survey galaxy that actually en- 
ters the 2dF the spectrograph is considerably less than the 
isophotal (or total) flux, as the 2dF instrument's optical fi- 
bres are only 1 arcsec in effective radius. The fibre magnitude 
of a galaxy of total magnitude m is given by 



where F g (< 9) is given in equation ([lB]). The z-dependence 
of mf is shown in in Fig. ^| for both an M* elliptical and an 
M* spiral galaxy. The higher central surface brightness of the 
elliptical results in a much higher fraction of its flux entering 
the 2dF fibre (dot-dashed line), but the PSF dominates the 
observed surface brightness of both types by z ~ 0.5, and so 
there is little difference in mi . On average mf — m' ~ 2 over 



most of the 2dF survey's redshift range (Section 4.5), which 
increases the survey's potential sensitivity to lensed quasars 
by ~ 2 mag. However this is somewhat offset by an analo- 
gous reduction in the fibre flux of any lensed quasars be hind 
the survey galaxies; this is discussed further in Section 6.1. 



mf 



m + 2.5 log 



F g (< t 



F g (< oo) 



(17) 



4.5 Redshift distribution 

The redshift distribution of the galaxies in the 2dF GRS is 
determined only by the (assumed) intrinsic galaxy popula- 
tion and the isophotal magnitude limit. In particular, is inde- 
pendent of the morphological selection criteria of the APM 
survey. The redshift distribution of a generic, magnitude- 
limited galaxy survey was discussed in Mortlock & Webster 
(2000b), and the results presented here were obtained by the 
same method, but with isophotal magnitudes m! replacing 
total magnitudes m throughout. Note that this incurs an 
additional computational overhead, as the conversion from 
isophotal magnitude (and redshift) to absolute magnitude 
requires an iterative solution. The resultant redshift distri- 
butions, dpg/dz, are shown as solid lines in Fig. |B| the dis- 
tributions that would be obtained if total magnitudes were 
used are shown as dashed lines. In both cases the narrower 
distribution was generated using the A'-corrections of Kin- 
ney et al. (1996); the broader distributions were obtained 
by assuming that the 7^-corrections are cancelled by pas- 
sive luminosity evolution of galaxies. The latter provides 
a reasonably simple model of the local galaxy population 
that reproduces the observed counts to m ~ 23, and also 
matches the observed redshift distribution of the 2dF galax- 
ies (e.g. Folkes et al. 1999) . The use of isophotal magnitudes 
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Figure 4. The total magnitude, m (solid line), APM isophotal magnitude, ml (dashed line), and the fibre magnitude, nif (dot-dashed 
line), of an M, elliptical galaxy (a) and an M* spiral galaxy (b), shown as a function of redshift, z. It is assumed that the passive 
luminosity evolution evolution of the galaxies cancels out their i\-corrections. The horizontal dashed line is the isophotal magnitude 
limit of the 2dF survey, ml,. = 19.5. 
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Figure 5. The expected redshift distribution, dp g /dz, of the 
galaxies in the 2dF GRS. The solid lines show the distributions 
generated if the correct isophotal magnitude limit of m[ im = 19-5 
is used; the dashed lines show the distributions implied by a total 
magnitude limit of mi; m = 19.5. In both cases the broader dis- 
tribution assumes that the i\-corrections are exactly cancelled 
by passive luminosity evolution of the galaxies; the narrower 
distribution assumes no luminosity evolution, and uses the K- 
corrections measured by Kinney et al. (1996). 



in place of total magnitudes has almost exactly the same 
effect as omitting the luminosity evolution, relative to the 
default model used in Mortlock & Webster (2000b). The 
2dF survey thus has the same depth as a GRS with a to- 
tal magnitude limit of m ~ 19 (as opposed to 19.5). Using 
the generic results of Mortlock & Webster (2000b), this re- 
duction in depth then implies a one third reduction in the 
number of lenses in the 2dF survey. 



5 MERGED OBJECTS 

Images separated by ^ 5 arcsec are treated as single objects 
by the APM software, and thus Maddox et al. (1990a) in- 
troduced the ^apm-Mapm plot to distinguish merged images 
from isolated stars and galaxies. The ide as that underpin 
this method were explained in Section 2.3; here the popula- 
tions of mergers shown in Figs. |l] and g are defined. 

A close pair of images of specified types are defined by 
their angular separation, A8, and their total magnitudes^*], 
mi and m,2. For chance alignments, the magnitudes are un- 
correlated, and so drawn at random from the respective pop- 
ulations, and the distribution of angular separation is given 
by Ap/A/S.9 oc AO between and A# max ~ 5 arcsec. The 
more complex separation distributions of binary stars and 
physically-associated galaxies are ignored here, as it is only 
the region of phase space inhabited by merged images that 
is of interest. 

Merged image pairs are shown in Figs. |l] and |^ as open 
symbols. Star-star merges are shown as open circles; star- 
galaxy mergers as open squares and galaxy-galaxy merg- 
ers as open triangles. All three populations have V'apm 
1000, and lie a significant distant from the stellar locus in 
Fig. pj. Whilst they also tend to be slightly below the galac- 
tic population, they would be treated as extended objects in 
the APM survey if not for the use of the &apm~/^apm plot 
(Fig. ^). Not only are the mergers are quite well separated 
from the monolithic objects (i.e. single galaxies and stars 
with /capm — 0), but the stellar images (either isolated stars 
or stellar mergers) are quite distinct from those with strong 
extended components. This separation allowed Maddox et 
al. (1990a) to reject most stellar mergers from the galaxy 



** Despite the fact that isophotal magnitudes are the 'natural' 
choice for simulations of the APM survey, the separate isophotal 
magnitudes of two merged objects are somewhat ill-defined. The 
total flux of an image pair is simply the sum of the two compo- 
nents' total fluxes, but their combined isophotal flux cannot be 
determined from their individual isophotal fluxes. 



© 2000 RAS, MNRAS 000, 



8 D. J. Mortlock and R. L. Webster 



sample, whilst retaining most of the galactic mergers. Un- 
fortunately, many gravitational lenses would also have been 
rejected, as discussed below. 



6 LENSED QUASARS 

The determination of the number of lenses that enter the 
2dF GRS splits neatly into two parts: a reasonably standard 
calculation of the lens population; and the application of the 
survey selection effects. 

The lens calculation uses the same models and meth- 
ods discussed in Mortlock & Webster (2000b): the deflectors 
are assumed to be a non-evolving population of isothermal 
spheres (with the possibility of both a core radius and ellip- 
ticity accounted for in Section ff^]); the quasar population 
is chosen to match the number counts of Boyle, Shanks & 
Peterson (1988) and the redshift distribution of the Large 
Bright Quasar Survey (Hewett, Foltz & Chaffee 1995); and 
conventional methods are used to evaluate the lensing prob- 
ability (e.g. Kochanek 1996). 

More complex, and more interesting, are the effects of 
the various APM selection effects on the lens statistics. To 
enter the 2dF GRS, a lensed quasar must satisfy the follow- 
ing conditions (c.f. Mortlock & Webster 2000b): 

(i) The composite object defined by the quasar images 
and lens galaxy must have an isophotal magnitude m' < 
m( im = 19.5. This is not a simple function of galaxy and 
quasar magnitudes, but must be computed numerically for 
each lens. Fortunately the quasar images are within the 
bright (i.e. /i <C fium) 'core' of the galaxy in many cases, 
and their total - rather than isophotal - flux can be added 
to the isophotal flux of the galaxy. 

(ii) The quasar must be bright enough to for its spec- 
troscopic signature to be present in the composite spec- 
trum obtained during the 2dF GRS. Following Kochanek 
(1992), the quasar is assumed to be detectable if the fi- 
bre magnitudes of the quasar images and galaxy satisfy 
mf,q — m f,g 5: Anif.qg ~ 2. The explicit cal cula tion of the 



fibre magnitudes is discussed in Sections 4.4 and 6.1 



(iii) The composite object must be classed as a galaxy 
according to the APM star-galaxy separation scheme. This 
condition is taken to be satisfied if m g — m a < Am gq , where 
the value of Am™ is estimated in Section 6.2. Note that m„ 



is the total, lensed magnitude of all the quasar images. 

The small aperture size tends to result in a greater sensitiv- 
ity to the more centralised quasar images (Section 3.1), but 
the use of isophotal magnitudes decreases the effective depth 
of the survey, and the rejection of multiple point sources also 
tends to reduce the number of lenses (Section 3.2). The rela- 



tive importance of these competing effects cannot be judged 
a priori, but must be determined numerically, as discussed 
in Section |^ 



6.1 Fibre magnitudes 

Fig. ^ shows that only a small fraction of the light from most 
galaxies enters the small 2dF optical fibres. However, much 
of the light from a lensed quasar can also fall outside the 
fibre, both due to seeing effects and the fact that the quasar 



images are generally somewhat offset from the centre of the 
fibre. 

The flux from a lensed quasar which enters the fibre 
is calculated by integrating the surface brightness of each 
quasar image over the fibre, and summing the resultant fibre 
fluxes. Thus the fibre magnitude, mf iq , of a lensed quasar 
with source position /3 is 



mf. 



/ Ei At [QiQ9)] T[0iQ9)] 

2 - 51os 1 I>W)1 



(18) 



where T(9) is the fraction of the light from an image at 9 
which enters the fibre. The image positions, 6i, and mag- 
nifications, are functions only of source position, /3, and 
are regarded as inputs to this calculation. Performing the 
two-dimensional integrals over the aperture can be time- 
consuming, but, for simple PSFs [including equation ([?])], 
the problem can be reduced to a one-dimensional integral 
(Mortlock 1999). 

The dependence of T on both 9i and 9 S is shown in 
Fig. [| The solid line in each plot shows the step function 
obtained if the PSF is ignored. The other four curves show 
that an appreciable portion of the flux from images up to 
~ 2 arcsec from the centre of the fibre is registered in the 
worst seeing to be expected at the AAT. This effect is more 
important for PSFs with more pronounced wings, as can be 
seen by comparing the results for the 'default' Moffat (1969) 
profile (r) — 2.5) shown in Fig. ^| (a) with those for Gaussian 
seeing (77 — > 00) shown in (b). 

Irrespective of the details of the PSF, Fig. ^ illustrates 
how the positions of the quasar images can influence rrif. q , 
and thus the chance of a lens being detected in the GRS. 
Obviously the OiS are determined by /3, but the types of 
image configurations produced depend on the mass distr- 
bution of the deflector. Several representative examples of 
this are shown in Fig. (?], in which the surface brightness of a 
galaxy-quasar composite is shown for a number of plausible 
lens models. The galaxy is an M* elliptical at a redshift of 
0.1, and has m g = 17.9 and m' g — 18.5, ignoring the quasar 
light. In all panels the quasar is slightly mis-aligned with 
the galaxy, and has an unlensed magnitude of 22. The as- 
sumption of 1 arcsec seeing means that the quasar images 
are sufficiently smeared that they have only a gross influ- 
ence on the surface brightness contours, and often remain 
unresolved. 

The default mass model, a singular isothermal sphere 
(as used for most of the calculations discussed in Section ^), 
is shown in (a). The two images lie just outside the fibre, 
but the PSF is sufficiently broad that the fibre flux of the 
quasar is three times its unlensed flux. 

In (b) the galaxy is no longer singular - the density 
flattens off within r c ~ 100 pc of its centre, resulting in 
the formation of an extra image and generally higher mag- 
nifications. This results in an increased fibre flux, but, for 
reasonable core radii, the effect is only minimal - just 10 
per cent in the example shown. Note that there is no signif- 
icant change in the image separation for small r c - use of 
a self-consistent dynamical normalisation (Kochanek 1996; 
Mortlock & Webster 2000a) cancels out the reduced depth 
of the galaxy's central potential well. For very large cores, 
the image separation actually increases, and hardly any of 
the quasar's light would enter the fibre. As shown in Fig. ^ 
(a) this reduces the lensing probability to zero. The proba- 
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Figure 6. The fraction of the light, T, from a point-like image that enters a 2dF optical fibre as a function of 8, the angular separation 
between the centre of the image and the centre of the 1 arcsec radius fibre. A Moffat (1969) PSF is assumed, with rj = 2.5 in (a) and 
Tj — > oo (Gaussian seeing) in (b). In both panels T is shown for a several values of the seeing: 8 S = arcsec (solid line); 8 S = 1 arcsec 
(dashed line); 8 B = 2 arcsec (dot-dashed line); 8 B = 3 arcsec (dotted line); and 8 B = 4 arcsec (dot-dot-dot-dashed line). 



bility of lensing by a spiral galaxy has a similar dependence 
on the core radius, but for very different reasons. The dy- 
namical normalisation of spirals is taken from their rotation 
curves, which are essentially unaffected by the core struc- 
ture. Hence increasing r c results in reduced image separa- 
tions, and an even greater increase in the quasar's fibre flux. 
However the shrinking of the radial caustic dominates the 
statistics, again reducing the lensing probability to zero for 
large cores. 

More important than the core structure of the deflec- 
tor is its ellipticity. A mass profile of ellipticity of e can 
be modelled adequately by applying an external shear of 
magnitude 70 = e/3 (e.g. Kassiola & Kovner 1993), and 
Fig. (?) (c) shows a singular lens with 70 = 0.2, represent- 
ing a galaxy of slightly greater than average ellipticity (e.g. 
Keeton, Kochanek & Seljak 1997). The source position was 
chosen so that it lies inside the lens's tangential caustic, and 
the extra pair of images contribute significantly to the de- 
tectability of the lens. The total magnification of the source 
is ~ 16, and the quasar's fibre magnitude is 20.2. Extra 
magnification is a generic feature of quadrupole lenses (e.g. 
Kochanek & Blandford 1987; Keeton et al. 1997), and the 
inclusion of ellipticity in the calculation can increase the ex- 
pected number of lenses by up to 50 per cent [Fig. ^ (b)] . 

Finally, Fig. (?) (d) shows the combined effect of a core 
radius and finite ellipticity - relative to the singular model, 
the magnifications are slightly increased and there is an ex- 
tra image, which increases the quasar's fibre flux by a further 
20 per cent. It can also been seen that, in rough terms, the 
effects of core radius and deflector ellipticity decouple, and 
it is the shape of the lens that dominates the lens statistics. 

Note that in all four of the above examples the overall 
magnification is high enough to ensure the lensed quasar's 
fibre flux is greater than its total unlensed flux would have 
been. Nonetheless, for most lenses more than half the light 
of the lensed quasar images falls outside the fibre - typically 
m f,q~ m q — 1- This is a significant loss, but the effect is still 
much smaller than it is for galaxies, which have mf ig — m g ~ 



2 on average (Section 4.4). Thus the small diameter of the 
2dF fibres is favourable for spectroscopic lens searches. 



6.2 APM classifications 

The APM classification of a lensed quasar as stellar, merged 
or non-stellar is a function of its surface brightness, and 
hence sensitive to both the properties of the lens galaxy 
and the positions and magnifications of the quasar images. 
The number of independent variables that 'define' a given 
lens is too high to facilitate an analytic investigation, and 
so Monte Carlo simulations of a large number of source- 
deflector pairs were generated. For each lens, the three APM 
image classification parameters (V>apm, &apm an d /iapm, as 
described in Section ^^) were calculated, and the simulated 
lenses are plotted as filled circles in Figs, [l] and ^. 

As expected, almost all the lenses have i/'apm ^ 1000, 
indicating that their surface brightness profiles are signif- 
icantly different from those of stars. However the popu- 
lation of lenses is quite similar to the merger population 
discussed in Section ^. In particular, two-image lenses with 
high-redshift deflectors are almost indistinguishable from bi- 
nary stars. Hence the rejection of star-star mergers from the 
APM galaxy survey on the basis of their location in ^apm- 
A*apm space also results in many lenses being lost from the 
2dF GRS. In fact the only lenses that remain in signifi- 
cant numbers are those with bright (and hence nearby) lens 
galaxies. Almost all lenses with m g — m q £ — 1 are morpho- 
logically classified as galaxies, and those with m g — m q <; 1 
are flagged as stellar mergers. In the calculation that fol- 
lows, this is simplified to the requirement that all lenses 
m a < Am EQ = ± 1. 



7 RESULTS 

The aim of this calculation is essentially to find one value: 
N\, the number of lensed quasars expected in the 2dF GRS. 



© 2000 RAS, MNRAS 000, 



10 D. J. Mortlock and R. L. Webster 





Figure 7. Simulated images of a a, = 220 km s _1 elliptical galaxy (at z g = 0.1 and with m g = 17.9) lensing a m q = 22 quasar at 
2 q = 2. The galaxy has m' ~ 18.5 in each of the panels, which are labelled with the fibre magnitudes of the galaxy, mf g , and the quasar, 
mf q . In all cases the quasar is off axis by ~ 0.2 arcsec, and its images are shown by the black disks, the areas of which are proportional 
to their magnifications. The heavy circle is the 2dF fibre, and the contours are uniformly spaced in logarithmic surface brightness. In 
(a) the galaxy is a singular isothermal sphere; in (b) it has a core radius of r c = 100 pc; in (c) there is an external shear of 70 = 0.2 
(equivalent to an ellipticity of e ~ 0.6); and in (d) the galaxy has both a core radius of r c = 100 pc and an external shear of 70 = 0.2. 
In all cases the seeing is 9 B = 1 arcsec, and it is assumed that the night sky has zero surface brightness. 



This is given by integrating the lensing probability of indi- 
vidual galaxies over the entire galaxy population, subjec t to 
the selection effects described in Section ^. In Section |7.l| 
the dependence of Ni on the survey design is investigated; 
the variation with lens model is discussed in Section W/A. 



7.1 Survey design 

The generic calculation of lensing in redshift surveys pre- 
sented in Mortlock & Webster (2000b) implies that ~ 10 
lenses would be discovered in a survey of Mot = 2.5 x 10 5 
galaxies that is complete to a total magnitude limit of 
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m = 19.5. Implicit in this calculation were the assumptions 
that there are no surface brightness-related selection effects 
in the survey, and also that the apertures used to gather 
the survey spectra are larger than the typical angular scales 
of both the galaxies and lensed images. These assumptions 
are dropped here, and Fig. ^ shows how Ni varies with var- 
ious aspects of the survey design; the generic model can be 
reproduced by setting jj,u m — > oo and Of — > oo. In this lim- 
iting case, both the fibre magnitude, nit, and its isophotal 
magnitude, m' , of an object approach its total magnitude, 
m. 

Fig. |^ (a) shows how N\ increases with fj-um, the isopho- 
tal limit of the survey. Most galaxies with m ~ 19.5 have 
observed peak surface brightnesses close to 22 mag arcsec -2 , 
and so the properties of the redshift survey (and N\ in par- 
ticular) do not vary with /iiim if the isophotal limit is either 
much smaller than or much greater than this critical value. 
The increase of N\ with fin m is simply due to the increased 
mean reds hift of the survey galaxy population, as described 
in Section 4.5. However for /^ii m <; 30 the isophotal and total 
magnitudes of all the objects considered here are essentially 
the same, and N\ flattens off. 

Interestingly, in the limiting case of fium — > oo, Ni is 
twice the value obtained in Mortlock & Webster (2000b). 
This is due to the finite aperture size, as shown in Fig. ^ 
(b). Taking 9f — > oo halves N\, as the increase in the fibre 
flux of the survey galaxies is much greater than that of the 
more centralised quasar images. However, Ni, does not in- 
crease for 6f 5s 1 arcsec, as there are few lenses with image 
separations less than ~ 2 arcsec. This implies that the 2dF 
fibre size is close to optimal for a lens search: any smaller and 
the integration times required to obtain reasonable spectra 
of m ~ 19 galaxies would become prohibatively long; any 
larger and fewer lenses would be found. 

Also of interest is the variation of N\ with the seeing, 
which is illustrated by the five lines in each panel of Figs. ^ 
and ^. In most cases there is little difference between 9 B = 
arcsec and 9 S — 1 arcsec, primarily because 9 S ^ 9f in this 
range. However if 9 S J; Of the number of lenses decreases 
with the seeing, because so much of the flux from even cen- 
tred quasar images misses the fibres, as illustrated in Fig. ^| 
(The survey galaxies are extended intrinsically, and so the 
relative reduction in fibre flux is lower.) The one exception 
to this qualitative explanation is shown in Fig. ^| (b), in 
which the number of lenses decreases rapidly with decreas- 
ing fibre size when S = 0. This case can be understood 
geometrically, as the quasar images must lie inside the 2dF 
fibres to be registered at all (there being no 'spillage' from 
the PSF), whereas most low-redshift galaxies produce image 
separations slightly larger than the fibre diameter. 

Fig. ^ (c) shows an exponential increase in N\ with 
Amf iqg , which reflects the steepness of the quasar luminosity 
function. For Amf iqg ^ 2 the change in Amf iqg corresponds 
almost directly to the increase in the effective depth of the 
'lens survey'. It is also clear that Nj is more sensitive to 
Amf iqg than it is to any other model parameter. If a de- 
tailed analysis of the 2dF spectra revealed that Am.f iqg ~ 2 
was an underestimate of the quasars' 'spectral prominence', 
the expected number of lenses could easily be doubled. 

As shown in Fig. ^ (d), the Am gq dependence is sim- 
ilar, but less pronounced than the dependence on Am.f !qg 
discussed above. This is because a high value of Am gq only 



results in the 'lens survey' being more complete to a given 
depth, whereas increasing Amf iqg probes fainter magnitudes 
and thus larger numbers of quasars. Combined with the fact 
that A m Kq is reliably constrained to be close to (from 



Section 6.2), the uncertainties in the star-galaxy separation 
properties of the APM survey do not flow through to the 
calculation of Ni. 



7.2 Lens model 

The most significant uncertainties in the calculation of N\ 
are relat ed to the observational parameters discussed in See- 
but the results do vary with the lens model, and the 



7.1 



tion 

effects of a finite core radius and ellipticity were explored. 

Fig. (a) shows how N\ varies with canonical core ra- 
dius, r c ». The core radius is assumed to vary with velocity 
dispersion as r c — r c »(<7 /a,) 4 , and r c * is taken to be the 
same for each galaxy type. The fall-off of N\ wit h co re ra- 
dius is very pronounced, as explained in Section 3.1 - for 
large r c » the spirals' cross-sections become zero and the im- 
ages formed by Es and SOs completely miss the 2dF fibres. 
So the assumption of a singular lens model overestimates Ni 
by at most 30 per cent, given that r c » < 100 pc (Kochanek 
1996). 

More important is the ellipticity of the lens galaxies, 
which are typically inferred to be as high as e ~ 0.6 from 
lensing studies (e.g. Keeton et al. 1997). This was investi- 
gated by the application of an external shear of magnitude 
70 = e/3, with all galaxies assumed to have the same pro- 
jected shape - a more realistic distribution would simply 
smooth the curves shown in Fig. |^ (b). As discussed in Sec- 



tion 6.1 



and shown in Fig. Q (c), the extra images that can 
be produced by quadrupole lenses almost always lie inside 
the 2dF fibres, and so the greater magnification bias can 
significantly increase N\. Given that 70 <; 0.3 is unrealistic 
for the dark matter-dominated galaxy halos, the use of a 
circularly-symmetric lens model underestimates the number 
of lenses by at most 50 per cent. 

Overall, the uncertainty in the deflectors' mass distri- 
butions does not greatly affect the estimates of N\ as the 
two effects described above tend to cancel out. The most re- 
alistic model would have a small (5s 50 pc) core radius and 
an ellipticity of between 0.3 and 0.6 (i.e. a shear of between 
0.1 and 0.2) which would result in only slightly more lenses 
than the default singular model. The calculation using sin- 
gular isothermal spheres also provides a robust lower limit 
of Ni ~ 10. 



8 CONCLUSIONS 

Kochanek (1992) and Mortlock & Webster (2000b) showed 
that appreciable numbers of lensed quasars could be discov- 
ered in galaxy surveys, but a number of potentially impor- 
tant selection effects and biases were ignored in both cal- 
culations. These generic results have been refined here to 
obtain a more realistic estimate of the number of lenses, N\, 
expected in the 2dF GRS. 

In particular, the importance of a finite collection area 
for the redshift survey spectra and a finite surface bright- 
ness limit were investigated. The 2dF optical fibres are un- 
usually small (an angular radius of only 1 arcsec), which 
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Figure 8. The expected number of lenses, JVj, in a GRS of the same design as the 2dF survey, as a function of limiting surface brightness, 
IHim, in (a), fibre radius, 9{, in (b), Am,f qg in (c) and Am gq in (d). The vertical dashed lines indicate the values appropriate to the 2dF 
GRS: /iii m = 24.5, 9{ = 1 arcsec, Amf qg = 2 and Am,f qg = 0. The lens model used was the singular isothermal sphere with no external 
shear. The results are shown for several values of the atmospheric seeing: 9 B = arcsec (solid lines); 9 B = 1 arcsec (dashed lines); 9 B = 2 
arcsec (dot-dashed lines); 9 B = 3 arcsec (dotted lines); and 9 B = 4 arcsec (dot-dot-dot-dashed lines). 



approximately doubles the number of lenses, as so much of 
the light from the extended galaxies misses the fibres. Con- 
versely, the use of isophotal rather than total magnitudes 
reduces the effective depth of the survey (by ~ 0.5 mag in 
the case of 2dF) which reduces the depth of the 'lens survey' 
by a similar amount. 

The 2dF input catalogue was determined by the APM 
star-galaxy separation algorithm, and so this was applied to 
a simulated population of quasar lenses, with the result that 
most were classified as multiple stellar images. Nonetheless, 
any lens in which the deflector galaxy is at least as bright as 
the (magnified) quasar images should enter the APM galaxy 
survey and hence the 2dF GRS. Importantly, lenses with 
nearby deflectors (e.g. Q 2237+0305) are especially valuable, 
as discussed in Mortlock & Webster (2000b). 

Whilst each of the above selection effects can change N\ 
by a factor of several, they tend to cancel each other out. 
The generic estimate in Mortlock & Webster (2000b) im- 
plied that between ten (if only lenses with bright deflectors 
were included) and fifty (if lenses with high-redshift deflec- 
tors were not rejected from the survey) new quasar lenses 



could be discovered in the 2dF GRS. Clearly it the former 
figure which is most relevant, and thus the prior prediction 
of ~ 10 low-redshift lenses is confirmed by the more detailed 
simulations presented here. The random uncertainty in this 
value is mainly due to shot noise, but the variability of the 
atmospheric seeing at the AAT site is also an issue. The 
most important systematic error is the quality of the survey 
spectra (characterised by Amf iqg ), as it is not clear how ef- 
ficiently galaxy spectra can be searched for the presence of 
quasar emission features. 

Irrespective of the above uncertainties, it is clear that 
redshift surveys, and the 2dF GRS in particular, are poten- 
tial sources of large numbers of lensed quasars. The next 
step then is to implement a systematic search for quasars' 
emission lines amongst the 2dF galaxy spectra. A notable 
success in this field was the spectroscopic discovery by War- 
ren et al. (1996) of a lensed emission line galaxy. In a related 
project, Willis et al. (2000) have developed a more rigorous 
technique, using template-subtraction to search for lensed 
emission line galaxies in a 2dF sample of field ellipticals. 

Unfortunately these methods are probably not applica- 
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Figure 9. The expected number of lenses in the 2dF GRS, N\, as a function of the the canonical core radius of the lens galaxies, r c «, 
in (a) and the external shear, 70, in (b). The calculation is performed assuming Amf qg = 2 and Am gq = in all cases. The results are 
shown for several values of the atmospheric seeing: 8 S = arcsec (solid lines); 8 S = 1 arcsec (dashed lines); 8 S = 2 arcsec (dot-dashed 
lines); 8 B = 3 arcsec (dotted lines); and 9 B = 4 arcsec (dot-dot-dot-dashed lines). 



ble to the 2dF spectra as the target galaxies are not as ho- 
mogeneous as those in the above sample. The most obvious 
choice of search technique is principal components analysis 
(e.g. Murtagh & Hecht 1987), as it is already being used to 
analyse the 2dF galaxy spectra (Glazebrook, Offer & Deeley 
et al. 1999; Folkes et al. 1999). Further, an important part 
of the existing analysis software is the flagging of 'unusual' 
spectrap], and it is a subset of these that will become the 
lens candidates. Even if there are thousands of 'unusual' ob- 
jects, it would be feasible to search for the presence of quasar 
emission features by eye, which would probably yield several 
hundred genuine lens candidates, most of which would then 
be easily rejected on morphological grounds. The complete- 
ness of the final candidate sample may well be quite low, 
but its efficiency (i.e. the fraction of lens candidates which 
are real) should be very high relative to conventional lens 
surveys (Mortlock & Webster 2000b) . The discovery of even 
one lens in this manner would not only show the worth of 
such searches, but be valuable in its own right. 
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tt At the time of writing, about 10 5 survey spectra had been 
obtained, reduced, and processed. No lenses or lens candidates 
have been reported, but little effort has been put into the analysis 
of irregular objects thus far. It is not clear whether the pipeline 
reduction algorithms will be sufficient to generate a useful list of 
lens candidates; this should become more apparent as the survey 
nears completion. 
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